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help of the dimensionless quantities. The dimensionless Key words: Oscillatory flow, 
partial differential equations were reduced to ordinary Double-Diffusion, Heat Source, 
differential equations with the aid of perturbation parameters, ODE, Concentration, and Velocity 
and the ODEs have been solved analytically where the profile, Temperature profile, 


temperature, concentration, and velocity profiles were Convection, Lewis Number. 


obtained. In the course of scaling and solving the governing 
system of the equation, we obtained some pertinent 
parameters with which we carried out numerical simulation 
using Wolfram Mathematica by varying the parameters such 
as thermal Rayleigh number, solutal Rayleigh number, Lewis 
number, Darcy number, magnetic field parameter, Prandtl 
number, heat source parameter, and oscillatory frequency 
parameter within a bounded domain. In conclusion, it is seen 
in the research that varying the pertinent parameters has an 
impact on the various flow profiles. 


1. Introduction 


When a fluid saturated porous layer confined between two infinite parallel plates containing two substances 
with different rates of diffusion is heated and salted below, a threshold temperature is reached where the 
warmer fluid at the bottom expands and becomes lighter and moves to replace the colder fluid at the top, 
thus creating a convection current. These currents originate when a body force acts on the fluid due to 
density gradients. The force which induces this convection current is called buoyancy, and the phenomenon 
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is called double-diffusive or thermosolutal convection. Copious literature has been devoted to studying this 
phenomenon. Pioneer studies on the subject are credited to Lord Rayleigh in his trail-blazing experiment in 
1916. Chandrasekhar [1] advanced conditions leading to the onset of stability for both stationary and 
oscillatory convection for free-free, free-rigid, and rigid-rigid boundaries. To predict the conditions leading 
to the onset of instability, several methods have been applied to certain degrees of accuracy, which validates 
experimental results. Some of these methods include Galerkin, normal mode analysis, perturbation, and 
energy methods. Many of the industrial applications of double diffusion convection include 
electrochemistry, crystal growth, petroleum reservoirs, magma chambers, insulation engineering, coal 
combustors, nuclear reactors, food processing, reservoir modelling, thermal insulation [2-3]. 


Comprehensive studies have been done on double-diffusive convection and can be found in books by 
Ingham and Pope [4], Nieldand Bejan [5], Batchelor [6], Drazin & Reid [7], Charru [8], Schmid & 
Heningson [9], Horton and Rogers [10], Wooding [11], Vafai [12]. Hill [13] analysed double diffusive 
convection in a porous medium with a concentration-based internal heat source using the Darcy and 
Boussinesq approximation. In this study, it was assumed that the porous boundaries are insulated, and the 
magnetic field is assumed to be small. The equation of state is taken to be linear with temperature and 
concentration. When the thermal Rayleigh number is less or equal to the heat source parameter, no 
oscillatory convection occurs. Compared to comparing the linear and nonlinear thresholds using the linear 
stability analysis and energy method, the findings showed that irrespective of the magnitude of the heat 
source parameter, the threshold remains unchanged even when the solutal Rayleigh number is increased. 


Gaikwad and Dhanraj [14] have examined the onset of double-diffusive reaction-convection in an 
anisotropic porous layer with an internal heat source. In this study, the boundaries of the wall are subjected 
to chemical equilibrium and the effects of thermal Rayleigh number, Damkohler number, thermal 
anisotropic parameter, internal heat parameter, and mechanical anisotropy parameter are to destabilize the 
system, thereby hastening the onset of instability, while the system is stabilized in the presence of thermal 
anisotropy parameter for both the stationary and oscillatory modes. Similarly, the presence of Damkohler 
number has a destabilizing effect on the stationary mode, while for the oscillatory mode; it delays the onset 
of instability and stabilizes the system. 


Israel-Cookey et al. [15] have studied the simultaneous effects of a vertical magnetic field and a 
concentration-based internal heat source on the onset of double-diffusive convection in a horizontal porous 
medium using the linear stability analysis technique. The study revealed the Hartmann and solutal Rayleigh 
numbers delay the onset of instability for both stationary and oscillatory convection, while the Lewis number 
stabilizes the system for stationary and hastens instability only for the oscillatory mode. The Lewis number 
destabilizes the system for oscillatory convection, whereas heat source parameters lead to a destabilization 
of the system for both stationary and oscillatory convection. 


Hill [16] has investigated the thermal instability of double-diffusive convection in an inclined porous layer 
with a concentration-based internal heat source. It was shown in this study that stabilization of the system 
was achieved with an increase in the angle of inclination for the transverse rolls, the magnitude of the 
vertical solutal Rayleigh number and internal heat source parameter. 


Israel-Cookey et al. [17] examined the combined effects of Soret and magnetic field on thermosolutal 
convection in a porous medium with a concentration-based internal heat source. The study showed that the 
effect of magnetic field is to delay the onset of instability for both modes. Soret destabilizes the system and 
for ranges of the heat source parameter less than, the system is stable; whereas, for values greater than or 
equal to, the system becomes unstable for all values of Soret and magnetic field. 


Kumar et al. [18] studied the linear and nonlinear thermosolutal instabilities in an inclined porous layer with 
a concentration-based internal heat source. Subcritical instability regions were obtained when the linear and 
nonlinear thresholds for longitudinal and transverse rolls were compared. Instability of the system is 
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hastened when the thermal diffusivity is greater than the solute and internal heat source. However, an 
increase in the angle of inclination stabilizes the system, which delays the onset of instability. 


Ebiwareme and Israel-Cookey [19] investigated the onset of Darcy-Brinkmann convection in a thin, porous 
layer with concentration-based internal heating with magnetic field effects utilizing the normal mode 
analysis technique. The findings of this study revealed that stabilization of the system is achieved in the 
presence of Vadasz and Darcy numbers, while the heat source parameter hastened instability in the system 
for both stationary and oscillatory convection. 


Most recently, Odok et al. [20] have considered the onset of magneto-convection in a rotating Darcy- 
Brinkmann porous layer heated from below with a temperature-dependent heat source using linear stability 
analysis for free-free boundaries. The effects of rotation, magnetic field, heat source, Prandtl number, and 
ratio of viscosities for both stationary and oscillatory convection are analysed and presented graphically. The 
study shows the rotation, magnetic field, and ratio of viscosities delay the onset of both stationary and 
oscillatory convection. The internal heat source parameter destabilizes the system while the onset of 
instability is slowed with the presence of Prandtl number. 


Our motivation in this present study is to extend the novelty of obtaining the velocity and temperature 
profiles for momentum, energy, and specie concentration equations by assuming oscillatory boundary 
conditions. The effects of the pertinent parameters on the different profiles are presented in tables and 
graphs. The organization of the study is as follows. Section 1 gives a detailed review of double-diffusive 
convection in a porous layer with a concentration-based internal heat source. The basic governing equations 
of the MHD flow subject to the prescribed oscillatory boundary condition are presented in Section Two. 
Section three gives the detailed mathematical analysis of the problem to give the different profiles for 
momentum, energy, and concentration profiles. The results of the obtained solution as presented graphically 
for different values of the model parameters are contained in section four. The discussion of the impact of 
the various parameters on the different profiles is given in section five. Concluding remarks are given in 
section Six. 


2 Mathematical Formulation 


We shall present the models governing the fluid flow through porous boundaries by following Issrael- 
Cookey et al.[17], which are 


2.1 Continuity Equation 
Vw=0 (1) 
2.2 Momentum Equation 


Ap, ow P“ * FB * C’\ ok 
Bi OW Fe +JxB-p(T*,C ) gk 2 
Ao oe k 


2.3 Energy Equation 
E E OT ‘ 
Aas twy )T = 07 a + O(C -C,) (3) 
2.4 Specie Concentration Equation 
É © (wy) =D, ad (4) 
ôt oy 


Boussinesq Approximation 
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p(T’.C’)=p,|1-B, T -T,)+Be(C’-C,) | (5) 
J=0,(Ē' +w xB’),V.J =0 (6) 
where A= G L ap = a = ia 

(ec,), (ee), (Pe), 


The boundary conditions are: 


w =0,T" -1,+(4-% Bc! -¢,+(S55) at y =h 


(7) 
w =0,T" =T, C =C, at y =0 


Assuming the fluid flow has no electrical conductivity, then E = 0,B = (0, B,,0) 
Following Israel-Cookey et al. (2017), the curl of the current and the magnetic field becomes 
i oe (w x B’) =-0 Bw (8) 


Using the following nondimensional variables 


faye ) Qr twa È w,o= a7 D) = Ra ES) 


h Ah? i E Ôc (9) 
ôc =C,—C,,6, =T, -T,, P= —(P" + ABV) E = A 
Ha, A 
Substituting equation (9) into equations (1)-(4), we have: 
V.w=0 (10) 
Da |ôw oP 2 
— =—-—+~VRad-M'w+Rs 11 
C A an 
06 180 
—=— —+NVRa 12 
ôt Le dy’ g es 
2 
p% Le K 
ôt Le ðy 
The subjected boundary conditions are: 
w=0,0=0,¢=0 at y=0 
J J (14) 
=0,0= m = = at y=1 


where the governing parameters becomes 
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Ra= Po& Êrkhô; Rs = Los Bckhoe M = ko By feu 22 
y Har Har Vu h? (15) 


2 
ya Pr= K Jas h esl py 


A Do Gy Da A h? 


3 Analytical Solution 


Since the flow is through an oscillatory porous boundary, It is assumed that the solution is in the following 
form: 


w=we”,0=He" b= pe” } (16) 


Differentiate equation (16) accordingly and substitute the result into equations (11)-(14), we have 


| Bee m -vraa — Ro =-P (17) 


Pro, 


2 
D — 8,0, =—NRaLeye ” (18) 
where p, = aLe 
8? 
st Bh =0 (19) 


where p, = ewLe 


The subjected boundary conditions are: 


m =0,4 =0,@ =0 at y=0 
J a (20) 
m= 0, 6, = = ed, = = ew at y =] 


Solving equation (19), we have 
p) = Ae) + Be) Q1) 


Solving equation (21) using the boundary condition in equation (20), we obtain the coefficients as: 


VJRae” JRae 


A=... 8 = —__ 22 
AAA e] a 
We substitute equation (22) into equation (21), so that it becomes: 
_|_ VRae*_| va) Rae” | (Wh) 23 
hl) (aM e] i a(i a] a 
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To get the temperature distribution, we shall solving equation (18), and the result obtained 
m J RaLey oo 
pz 


Solving for the constant coefficients in equation (24) using the corresponding boundary conditions in 
equation (20), which is 


B VRaLey oo 


(y) = 4e + Be”? (24) 


A, +B, = (25) 
' B, 
Ajev® + Be VP aj) NE ROT e” (26) 
2 B, 
Presenting equation (18) in matrix form, we have: 
1 1 \(A, _ VRaLey oo 
p 
= (27) 
oh Pale B VRa NRalLey ) n 
2 —— —--_—_—* |g 
2 p 
Let a, = ore and @, -( ži -— — j so that equation (27) reduces to: 
2 2 
1 1 \(A,) (a, 
= (28) 
eve e || B a, 
1 1 
Aye = (P -— gl] (28) 
o eve 
A = 7 pak E ae VP -a (29) 
a, e VP 
1 a, 
A, = = a, — eP (30) 
ovh: as 
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A, ae = snap = aE G1) 
2 
ard aa aa) 
Substitute equation (31) into equation (24), we obtain: 
ae” eos Bey Q, — a eP -JF P \ RaLey oo (32) 


0,(y)= (FF) ew + (e-e) e B. 


To obtain the fluid velocity profile, shall substitute equations (23) and (32) into equation (17), which is 


To obtain the fluid velocity, we substitute equations (23) — (32) into equation (17), we obtain: 


where all constants coefficients in equation (33) remain as mentioned above. 


Solving for the velocity, we have 


ae [as P 4. Boe Nh {Palen J; 
eae 
nee (34) 


ae [ae A) BE “H-5 
(Bae (Bee ae 
Pro, Pré 


But if we substitute the constants into equation (34), we have: 


( peo, a hw 7 JRa ae OA o A a, -aet oA» z Vv RaLey oot 
0 
(e-e) (eV — oF | B, 


(35) 


Rae“ (Ay) Rae“ (VR) 
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Solving for the perturbed velocity profile in equation (35), we have: 


P VRa ane mec) e Bay } a, -o,e e VA + VRaLey em 

(=) 2 (=) 2) (eM -e*) (e-e) 2B 

+M +M 
PeR Pro (36) 
Wy = 

Rs VRae” (VA) Rae” -( V>) 

(22) m>) 2d -e 0P) | (Ai) dm) 

Pré, 


Having obtained the perturbed solutions for concentration, temperature and velocity profiles, we shall 
substitute equations (24), (32), (35) and (36) into equation (16), which is 


_ Rae (Jy) n Rae (Jy) a9 
At) T) ° a(d] Á 
-JA _ ay enh 
6,(y)= Qe a, en PY y A, — Qe ov , VRaLey ow (38) 
(e-e) (e -e | hə 
P S VRa [as Poy + Be% Boy „ NRaLey e=] 
(Gaga [Bae p 
Tye Pro Pro, (39) 
= 
5 Rs (ad By) +B + aen 
Dao +M? 
Pro, 
Pe” | e” Ra ae VP =% |E ,| & -aef oe 4 VRaLey oe 
e)a) (pa ae) (ee) A 
Pré, Pré, (40) 
w(y,t)= 
| Rse™ VRae d By) | Rae“ A By) 
(22) n) 2d) A) afet) ~e) 
Pré, 
4 Results 


In this section, numerical computation was performed using computational software called Wolfram 
Mathematica, version 12. The simulation was done by varying some of the pertinent parameters involved in 
the analytical solutions. The parameters are derived from Israel-Cookey et al. [17], and are in ranges: 
O<M <5,0<y<1,0< Da<2,0< Pr<5, Rs=50 0< Le <3, 


© 2022, CAJMTCS | CENTRAL ASIAN STUDIES www.centralasianstudies.org ISSN: 2660-5309 | 32 


Copyright (c) 2022 Author (s). This is an open-access article distributed under the terms of Creative Commons 
Attribution License (CC BY).To view a copy of this license, visit https://creativecommons.org/licenses/by/4.0/ 


CENTRAL ASIAN JOURNAL OF MATHEMATICAL THEORY AND COMPUTER SCIENCES Vol: 03 Issue: 06 | Jun 2022 


10< Ra < 50. The graphical results are presented as follows: 


wiy,t) 


Fig 1 Velocity Profile with varying Magnetic Field and other values are 
Pr =0.72, Da = 0.5, Ra = 50, Rs = 50, Le = 0.5,t = 5,@ =3,y =5 


7.1588 = 10 


7.1587 x 10° 


wly,t) 


7.1586 x 10° 


7.1585 = 10° 


0.0 0.2 04 06 0.8 1.0 


y > 
Fig 2 Velocity Profile with varying Rayleigh Number and other values are 


Pr = 0.72, Da =0.5,M =0.2, Rs = 50, Le = 0.5,t = 5,@ =3,y =5 


7.1588 x 10° 


7.1587 x 10° 


w(y,t) 


7.1586 x 10° 


7.1585 x 10° 


Fig 3 Velocity Profile with varying Solutal Rayleigh Number and other values are 
Pr = 0.72, Da=0.5, Ra =50,M =0.2, Le =0.5,t=5,@=3,y =5 
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wily,t) 
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y—+ 
Fig 4 Velocity Profile with varying Heat Source and other values are 


Pr = 0.72, Da = 0.5, Ra = 50, Rs = 50, Le = 0.5,t =5,@=3,M =0.2 


Aly,t) 


0.0 0.2 04 06 08 1.0 
y > 
Fig 5 Temperature Profile with varying Rayleigh Number and other values are 
Le =0.5,t=5,@=3,y=5 
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0.0 0.2 0.4 0.6 0.8 1.0 


y > 
Fig 6 Temperature Profile with varying Lewis Number and other values are 
Ra = 50,t = 5,0 =3,y =5 
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Aly,t) 


Fig 7 Temperature Profile with varying Oscillatory Frequency and other values 
are Ra=50, Le =0.5,t=5,y =5 


ary, t) 


y > 
Fig 8 Temperature Profile with varying Heat Source and other values are 
Ra = 50, Le =0.5,t = 5,0 =3 


Hyt) 


Fig 9 Concentration Profile with varying Rayleigh Number and other values 
are Le =0.5,t=5,@ =3 
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iyt 


y > 
Fig 10 Concentration Profile with varying Lewis Number and other values are 
Ra=50,t=5,@=3 


Hyt) 


Fig 11 Concentration Profile with varying Oscillatory Frequency and other 
values are Ra = 50, Le =0.5,t = 5 


5 Discussion 


In Israel-Cookey et al. [17] investigated the stability of the system and studied the impact of the varying 
parameter values on temperature and concentration profiles. However, in this article we investigated the 
impact of the pertinent parameters on the flow profiles, including the velocity profile. In the analytical 
solution, it is seen that the velocity profile exhibits a linear relationship for all the pertinent parameters. After 
thorough simulation, we shall discuss the results as follows: 


Fig. 1 depicts the effect of the magnetic field on the fluid velocity. This result indicates that an increase in 
the magnetic field decreases the fluid velocity. This result typically shows the linear relationship of the 
parameters earlier mentioned. The Rayleigh number impact on the fluid velocity was investigated and the 
result is seen in Fig 2. This result is of the view that any increase in Rayleigh number increases the fluid 
velocity. The result of the impact of solutal Rayleigh number on the fluid velocity is shown in Fig 3. The 
figure shows that fluid velocity increases for an increase in solutal Rayleigh number. Figure 4 depicts the 
effect of a heat source on fluid velocity; this result indicates that fluid velocity increases as heat from a 
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source increases. The fluid involves the effect of the pertinent parameters on the fluid temperature; Fig 5 
depicts an increase in fluid temperature for an increase in Rayleigh number. Fig. 6 shows a rise in fluid 
temperature for every increase in Lewis number. This result showed that the temperature is zero when the 
boundary condition is zero, and the temperature grows to a maximum level around the walls due to wall 
thickness. Fig. 7 shows that the fluid temperature decreases for an increase in oscillatory frequency with 
different maximums of the profiles. The fluid temperature increase is seen in Fig 8. This figure shows that 
the fluid temperature increases for an increase in heat through the source. This research also includes the 
impact of the pertinent parameters on the concentration of the fluid. Fig 9 demonstrates that the 
concentration of particles increases for an increase in Rayleigh number from 10 to 50. Fig 10 showed that 
the concentration in the fluid decreases for an increase in Lewis number. Finally, it is seen in Fig 11 that the 
specie concentration in the fluid causes an increase in oscillatory frequency. 


6 Conclusion 


The exact solutions of the governing equations were obtained by reducing them to ODEs using the perturbed 
function. The numerical simulation was carried out using Wolfram Mathematica, version 12, where the 
results were presented graphically. We can conclude from the results as follows: 


1. The fluid velocity can be slowed down if we increase the magnetic field intensity. 


2. The fluid velocity is increased if the solutal Rayleigh number and thermal Rayleigh number are 
increased. 


When heat is passed through a source, the fluid velocity increases. 


4. The fluid temperature can be increased if we increase the thermal Rayleigh number, the Lewis number, 
and the source heat, while the temperature reduces for an increase in the oscillatory frequency. 


5. The specie concentration increases for an increase in thermal Rayleigh number, while the concentration 
reduces for an increase in Lewis number and oscillatory frequency. 
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